Probing the Dark Energy with Quasar Clustering 
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We show, through Monte Carlo simulations, that the Alcock-Paczyhski test, as applied to quasar 
clustering, is a powerful tool to probe the cosmological density and equation of state parameters, 
fl ra o, QxO and w. By taking into account the effect of peculiar velocities upon the correlation function 
we obtain, for the Two-Degree Field QSO Redshift Survey (2QZ), the predicted confidence contours 
for the cosmological constant (w — — 1) and spatially flat (f2 m o + QxO = 1) cases. It turns out that, 
for w = — 1, the test is especially sensitive to the difference fi m o — Oao , thus being ideal to combine 
with CMB results. We also find out that, for the flat case, it is competitive with future supernova 
and galaxy number count tests, besides being complementary to them. 
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Introduction. Recent investigations of type la super- 
novae (SNe la) suggest that the expansion of the Uni- 
verse is accelerating, driven by some kind of negative- 
pressure dark energy |Q,|). Independent evidence for 
the SNe la results is provided by observations of cosmic 
microwave background (CMB) anisotropies in combina- 
tion with constraints on the matter density parameter 
i^mo) Hj- The exact nature, however, of this dark en- 
ergy is not well understood at present. Vacuum energy 
or a cosmological constant (A) is the simplest explana- 
tion, but attractive alternatives like a dynamical scalar 
field (quintessence) Q have also been explored in the lit- 
erature. An important task nowadays in cosmology is 
thus to find new methods that could directly quantify 
the amount of dark energy present in the Universe as 
well as determine its equation of state and time depen- 
dence. New methods may constrain different regions of 
the parameter space and are usually subject to differ- 
ent systematic errors, and they are therefore crucial to 
cross-check (or complement) the SNe results. 

The test we focus on here is the one suggested by Al- 
cock and Paczyhski (hereafter AP) ||, which has at- 
tracted a lot of attention during the last years [(3| [lOf . 
In particular, Popowski et al. Q| (hereafter PWRO) ex- 
tended a calculation by Phillips JlJ] of the geometrical 
distortion of the QSO correlation function. They sug- 
gested a simple Monte Carlo experiment to see what con- 
straints should be expected from the 2dF QSO Redshift 
Survey (2QZ) and the Sloan Digital Sky Survey (SDSS). 
However, they did not estimate the probability density in 
the parameter space and, as a consequence, they could 
not notice that the test is in fact very sensitive to the 
difference f2 r „o — 17 ao- Further, they did not take into 
account the effect of peculiar velocities, although they 
discussed its role arguing that it would not overwhelm 
the geometric signal. 

Our aim, in this Letter, is to show the feasibility of 
redshift distortion (geometric + peculiar velocity) mea- 
surements to constrain cosmological parameters, by ex- 



tending the PWRO Monte Carlo experiments and obtain- 
ing confidence regions in the (r2 TO o,f2Ao) an d (fimOiw) 
planes. We compare the expected constraints from the 
AP test, when applied to the 2QZ survey, with those 
obtained by other methods. We include a general dark 
energy component with equation of state P x = w p x , with 
w constant. Our analysis can be generalized to dynami- 
cal scalar field cosmologies as well as to any model with 
redshift dependent equation of state. Since most quasars 
have redshift z < 2 we expect the test to be useful in 
the determination of a possible redshift dependence of 
the equation of state. We explicitly take into account 
the effect of large-scale coherent peculiar velocities. Our 
calculations are based on the measured 2QZ distribution 
function and we consider best fit values for the ampli- 
tude (rrj) and exponent (7) of the correlation function 
as obtained by Croom et al. |l5| ]. In this work, we only 
consider the 2QZ survey although the results can easily 
be generalized to SDSS. 

Alcock-Paczyriski test and quasar clustering. We as- 
sume that the geometry is described by the standard 
Robertson- Walker metric. By a straightforward calcu- 
lation for null geodesies, we obtain the radial coordinate 
R as a function of z: 

( sinh(Vn^I(z))/(H y/Tho), ^fco > 0, 
a B : = g(z) := < I(z), Ofcp = 0, 

[ sm{^h^i(z))/(H ^n^), n kQ < 0, 

(1) 

where ag is the present scale factor, I(z) := 
f z ,_ Q [Ho/H(z')]dz', and the Hubble parameter is given 
by H(z) = H Q p mQ {\ + zf + 0,0(1 + z) 3 ( 1+ -) + fifco(l + 
z) 2 } 1 / 2 . 

Given two close point sources (e.g., quasars), with co- 
ordinates (z, 9, (f>) and (z+dz, 6+d9, <f)+d<j)), directly read 
off a catalogue, the real-space infinitesimal comoving dis- 
tance between them can be decomposed, in the distant 
observer approximation we adopt, into contributions par- 
allel and perpendicular to the line of sight, r± :— g{z)da, 



r\\ := dz/H(z), such that r 2 



Here, da is the (z + dz, s + ds, fi + dfj,), 



small angle between the lines of sight. 

The gist of the AP test relies then on the fact that, 
if we observe an intrinsically spherical system (r| = 
rj_), it will appear distorted, in redshift space, accord- 
ing to the generic formula ri/rii = j(z)sx/su, where 
the anisotropy or distortion function j(z) is defined by 
j(z) :— g(z)H(z)/z. Here we have assumed a Euclidean 
geometry for redshift space, that is, sh := dz, s± := zda, 
and s 2 = s| + s 2 ± . 

Observations [|13| suggest that, on scales ~ 1 — 
40/i _1 Mpc, the real space correlation function for 
quasars is reasonably well fitted by a power law, 
£(r) = (t/vq) 7 , which leads, in redshift space, 
to an anisotropic correlation function, £(s, /i) = 

[s/so(2)]~ 7 [/i 2 +i 2 (^)(l-M 2 )] _7/2 , where /j := S] \/s 
and So (z) :— tqH(z). 

Peculiar velocities also induce distortions in the corre- 
lation function which can be confused with those arising 
from the cosmological geometric effect. It is important 
to take them into account when comparing theory with 
observations. For the (z, s) range we will consider, the 
influence of small-scale velocity dispersions is likely to be 
weak [|ll| and we neglect it in our analysis. The most 
relevant effect to be considered is due to large-scale co- 
herent flows [|14| . The linear theory correlation function 
is given by |1|§ 
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£(r), (2) 



where the Pi(n) are Legendre polynomials, and /1 := 
r\\/r. As usual, (3{z) := f(z)/b(z), /(z) := 
— d\nD/d\n(l + z) is the linear growth rate, and we 
adopt the following dependence for the bias parameter, 
b(z) = 1 + [D(z = 0)/D(z)} m (6 - 1) . If m = 1, we 



have Fry's number-conserving bias model |16|. The case 
m = corresponds to a constant bias, and we also use 
771 ~ 1.7 in our computations, which seems to be more in 
accordance with an observed nonevolving clustering p3[ . 
For models where the dark energy is a cosmological con- 
stant (w = —1), we use the Heath solution for the grow- 
ing mode 0, D(z) = |^ m0 iJ(z) /~ (1 + x)/H{xfdx 7 
and the following approximation for the growth 
rate @, /(z) ~ f>l /7 (z) + ^fi A (z) [1 + fi m (z)/2] . 
For flat models, Silveira and Waga jl9[ obtained 
an exact solution for the growing mode, D(z) = 



2F1 



1 



1 



1— n„ 



2 - „»' nm /(l + ^) 3t0 J/(l + ^where 
2-F1 [a, 6, c; x] is the hypergeometric function. The growth 
rate can also be expressed in terms of hypergeometric 
functions. 

Following PWRO, we obtain, for the number of pairs 
expected in an infinitely small bin within (z, s, /1) and 



dN, 



pairs 



2ir /180 N Q F{z) 
A\ irz 



[1 + £,l(s, m)] s 2 dzdsd/j,. 
(3) 



Here A is the area (in deg 2 ) of the survey, Nq is the total 
number of sources (quasars) in the survey, and F(z) is 
the normalized distribution function. 

Croom et al. [Q~3[ , assuming an Einstein-de Sitter Uni- 
verse (f2 m = l,f2A = 0)i showed that the quasar clus- 
tering amplitude r$ appears to vary very little over the 
entire redshift range of the 2QZ survey. They found 
ro ~ 4 ft," 1 Mpc as their best fit, which remains nearly 
constant in comoving coordinate. Therefore, we have 
so\ Ed s(z) := H(z)\ E dsr = (4/3000)(l + z) 3 / 2 . Following 
again PWRO, we use the fact that the total number of 
correlated pairs, N pa i rs , in the survey is model indepen- 
dent to scale sq(z) to other cosmologies. It is straight- 



forward to show that N„ 



ex Sq(z)/j (z), and we use 



i2/3 



s Q (z) = (4/3000)(l + z) 3 / 2 l3(z)/j\ E ds(z)Y /0 as a fiducial 
redshift-space correlation length for our simulations. 

A particular model predicts a number of pairs Ai in 
each bin of a (z, s, ^1) space. In a real (or simulated) sit- 
uation the data consist of A^ pairs in i bins. PWRO 
showed that for typical surveys, such as SDSS and 2QZ, 
we are bound to be in the "sparse regime" or "Poisson 
limit" . In this case we may treat bins in (z, s, /i) space as 
independent and the probability of detecting JVj pairs in 
bin i, when At are expected, is P(A l | J 4 l ) = e~ Ai Af i /Ni\ 
Since the bins are independent, the likelihood C of ob- 
taining the data given the model is simply the prod- 
uct, C = ]XP(iV;|Aj). For a typical 2QZ simulation, 
we assumed: (i) the completed survey will comprise 
Nq = 26000 quasars in a total area A = 750 deg 2 ; (ii) 
the Einstein-de Sitter fiducial correlation function has 
ro = Ah^ 1 Mpc and 7 = 1.6; (iii) the bias model is de- 
termined by bo = 1.5 and 777 = 1. The linear binning we 
chose covered the ranges: 0.4 < z < 2.6, 2 < s/s (z) < 7, 
and < /-J. < 1, with 16 bins in z, 25 in s/sq(z), and 5 
in ^t, making up a total of 2000 bins. The maximiza- 
tion of the likelihood was carried out with MINUIT |2Ql 
and cross-checked with MATHEMATICA. The probability 
density function was built via a Gaussian kernel density 
estimate, from typically 1000 runs for each "true" model. 

Results and discussion. In Figure [IL we show the pre- 
dicted AP likelihood contours in the (O m o, f2Ao)-pl a ne 
for the 2QZ survey (solid lines), in the case w = — 1, in 
a universe with arbitrary spatial curvature. The scat- 
tered points represent maximum likelihood best fit val- 
ues for Q m o and f^Ao- The assumed "true" values are 
(Q m0 = 0.3, fi A0 = 0) and (0.28,0.72), for the top and 
bottom panels, respectively. In the top panel the dis- 
played curve corresponds to the predicted 2a likelihood 
contour. In the bottom panel the predicted la contour 
(dashed line) for one year of SNAP data |2l]] is displayed, 



together with the predicted lex AP contour. For the 
SNAP contour, it is assumed that the intercept M. is ex- 
actly known. To have some ground of comparison with 
current SNe la observations, in the same panel, we also 
plot (dotted lines) the Supernova Cosmology Project ||] 
la contour (fit C). As expected, in both cases, the test 
recovers nicely the "true" values. We stress out that the 
test is very sensitive to the difference f2 m o — SIao- From 
the bottom panel we note that the sensitivity to this dif- 
ference is comparable to that expected from SNAP, of 
the order ±0.01. Comparatively, however, the test has 
a larger uncertainty in the determination of fi m o + ^A0j 
of the order ±0.21. The degeneracy in fi m o + Qao may 
be broken if we combine the estimated results for the AP 
test with, for instance, those from CMB anisotropy mea- 
surements, whose contour lines are orthogonal to those 
exhibited in the panels (22J . 




FIG. 1. Simulated models at fixed w — — 1 and correspond- 
ing predicted AP confidence contours (solid lines) . In the top 
panel we show the predicted 2a likelihood contour assuming a 
"true" model (f2 m o = 0.3, Qao = 0). In the bottom panel the 
predicted la contour (dashed line) for one year of SNAP data 
[Ell is displayed, together with the predicted la AP contour. 
For both tests we consider n m o = 0.28 and ^ao = 0.72; also 
displayed is a la confidence contour obtained by the Super- 
nova Cosmology Project (dotted lines; pi). 

In order to estimate the consequences of neglecting 
the effect of linear peculiar velocities, in the top panel 
of Figure @, we included them in the calculation of the 
Ai values but neglected them in the computation of the 
maximum likelihood; in this panel, we assume £l m o = 0.3 
and f^Ao — as "true" values. Notice that the point with 
the "true" tt m o and f^Ao values is outside the 2a contour. 
It is clear, therefore, the necessity of taking this effect in 
consideration when analyzing real data. 

To illustrate that the AP test is in fact more sen- 
sitive to the mean amplitude of the bias rather than 



to its exact redshift dependence, we plot, in the bot- 
tom panel of Figure @, the 2a contour line, assuming 
as "true" values f2 m o = 0.3 and f^Ao = 0.7. For this 
panel, the "true" Ai values were generated assuming 
bo = 1.45 and m = 1.68. However, for the simula- 
tions, we considered a constant bias (m = 0), such that 
bo,sim '■= J ™" x F (z)btrue{z)dz — 2.46. We remark that 
the contour is slightly enlarged and shifted in the direc- 
tion of the "ellipsis" major axis. However, the uncer- 
tainty in Q m o — ^ao is practically unaltered, confirming 
the strength of the test p3| . We did the same analysis 
assuming Q m o = 1 and f^Ao = and obtained similar 
results. 





FIG. 2. Simulated models at fixed w — — 1 and correspond- 
ing 2a predicted AP confidence contour; in both panels, the 
"true" model is indicated by a solid dot. Top panel: The 
"true" model, (0.3,0), takes into account the effect of pecu- 
liar velocities, but the simulated ones do not. Notice that 
the "true" model does not fall into the 2a confidence region. 
Bottom panel: The "true" model, (0.3, 0.7), uses a redshift de- 
pendent bias function with bo = 1-45 and m = 1.68, whereas 
the simulated ones use a constant bias equal to 2.46. 

In Figure pi we show the predicted AP likelihood con- 
tours in the (O m o,if)-plane for the 2QZ survey (solid 
lines) for flat models (flko = 0). The "true" values are 
(n mQ = 0.28, w = -1) and (fi ro0 = 0.3, w = -0.7) for the 
top and bottom panels, respectively. In the top panel, 
we show, besides the AP contour, the predicted contour 
for one year of SNAP data (dashed line; ]2l|j ), both at la 
level. For the SNAP contour, the intercept M. is assumed 
to be exactly known. Notice that the contours are some- 
what complementary and are similar in strength. In the 
bottom panel, we compare the predicted 95% confidence 
contour of the AP test with the same confidence contour 
for the number count test as expected from the DEEP 
redshift survey (dashed line; |24|]). Again the contours 
are complementary, but the uncertainties on f2 m o and w 
for the AP test are quite smaller. 
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FIG. 3. Simulated flat models and corresponding predicted 
AP confidence contours (solid lines). The top panel is from 
a "true" model (Q m o — 0.28, w — — 1), and displays the 
predicted confidence contours for the AP test and the SNAP 
mission (dashed line; pi|), both at la level. The bottom 
panel is from a "true" model (f2 m o = 0.3, w = —0.7), and 
displays the predicted confidence contours for the AP test 
and the DEEP survey (dashed line; 0), both at the 95% 
level. 

In summary, we have shown that the Alcock-Paczyhski 
test applied to the 2dF quasar survey (2QZ) is a potent 
tool for measuring cosmological parameters. We stress 
out that the test is especially sensitive to fi m o — f^AO- 
We have established that the expected confidence con- 
tours are in general complementary to those obtained by 
other methods and we again emphasize the importance 
of combining them to constrain even more the parame- 
ter space. We have also revealed that, for fiat models, 
the estimated constraints are similar in strength to those 
from SNAP with the advantage that the 2QZ survey will 
soon be completed. 

Of course our analysis can be improved in several as- 
pects. For instance, for the fiducial Einstein-de Sitter 
model, we have assumed that 7 and r^ do not depend 
on redshift. In fact, observations |l3f seem to support 
these assumptions, but further investigations are neces- 
sary. Further, in the simulations, for Figure FJ and Figure 
H, we have assumed that the parameters tq, 7, 60 and 
m are known exactly, that is, they are the same as the 
"true" input ones. Marginalization over these parameters 
is expected to increase the size of the contours. However, 
preliminary results where the errors in r$ and 7 are taken 
into account (supposed Gaussian), show that the confi- 
dence contours are not appreciably altered. At present, 
the quasar clustering bias is not completely well under- 
stood. Theoretical as well as observational progress in 
its determination will certainly improve the real capacity 
of the test. However, confirming previous investigations 
p3[, we have found that the test is, in fact, more sen- 



sitive to the mean amplitude of the bias rather than to 
its exact redshift dependence. A more extensive report 
of this work and further investigations will be published 
elsewhere. 
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